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The 'plate tectonics' is an observed fact and most models of earthquake incorporate that through 
the frictional dynamics (stick-slip) of two surfaces where one surface moves over the other. These 
models are more or less successful to reproduce the well known Gutenberg-Richter type power law 
in the (released) energy distribution of earthquakes. During sticking period, the elastic energy gets 
stored at the contact area of the surfaces and is released when a slip occurs. Therefore, the extent of 
the contact area between two surfaces plays an important role in the earthquake dynamics and the 
power law in energy distribution might imply a similar law for the contact area distribution. Since, 
fractured surfaces are fractals and tectonic plate-earth's crust interface can be considered to have 
fractal nature, we study here the contact area distribution between two fractal surfaces. We consider 
the overlap set (m) of two self-similar fractals, characterised by the same fractal dimensions (d/), 
and look for their distribution P(m). We have studied numerically the specific cases of both regular 
and random Cantor sets (in the embedding dimension d — 1), gaskets and percolation fractals (in 
d = 2). We find that in all the cases the distributions show an universal finite size (L) scaling 
behavior P'(m') = L a P(m, L); m' = mL~ a , where a = 2(df — d). The P(m), and consequently the 
scaled distribution P'(m'), have got a power law decay with m (with decay exponent equal to d) 
for both regular and random Cantor sets and also for gaskets. For percolation clusters, P(m) (and 
hence P'(m')) have a Gaussian variation with m. 



I. INTRODUCTION 

Enormous efforts have been made by geologists 
and physicists to understand the earthquake 
phenomena since several decades. Its dynam- 
ics is still a challenging problem. Though 
model studies of earthquake have been going 
on several decades now, there is no consensus 
regarding one single model. 'Plate-tectonics' 
is an important observation in this context 
by the geologists and the Gutenberg-Richter 
power law is an unique characterisation of this 
dynamics. The earth's solid outer crust (about 
20 km thick) rests on a tectonic shell. This 
tectonic shell is divided into a number (about 
12) of mobile plates, having relative velocities 
of the order of few centimeters per year. This 
motion of the plates arises due to the powerful 
convective flow of the earth's mantle, at the 
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inner core of earth. Due to the surface rough- 
ness of both the earth's crust and the plates, 
solid-solid frictional forces arise and this helps 
sticking the earth's crust to the plates. This 
sticking develops elastic strains and the strain 
energy gradually increases because of the uni- 
form motion of the tectonic plates. There is 
a competition between the sticking frictional 
force and the restoring elastic force (stress). 
When the accumulated stress exceeds the fric- 
tional force, a slip (earth quake) occurs and 
it releases the stored additional elastic en- 
ergy. Gutenberg and Richter (1954) pj, by 
analysing the records, observed a power-law 
distribution of the elastic energy released dur- 
ing earthquakes given as N(e) ~ e~ a , where 
N(e) is the number of earthquakes releasing 
energy greater than or equal to e and a is 
the power exponent. The observed value of a 
ranges between 0.7 and 1.0. Several hypothe- 
sis and model systems have been proposed to 
investigate the nature of the earthquake phe- 
nomena The main intention of the model 
studies is to capture the above Gutenberg- 
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Richter type power law for the frequency dis- 
tribution of failures (quakes) in the failure dy- 
namics of the models. Among these, most of 
the models 0,0 incorporate the stick-slip pro- 
cess and the roughness Q of the surfaces in- 
volved as important features. Recently some 
models 0, Q] have included the fractal nature 
of both the earth's crust and the plate involved 
in stick-slip process and initiated new mod- 
ellings of earthquake dynamics. As the elas- 
tic energy can be stored at the contact area of 
the surfaces only, the contact area distribution 
should have much importance. Also, in frac- 
tal physics the statistics of two fractal overlap 
is still missing although this may be useful to 
study the interface properties in many physi- 
cal situations. 

In this report we have studied the contact 
area distributions and their scaling properties 
for different sets of fractal surfaces, using com- 
puter simulation techniques. 



II. SIMULATION STUDIES OF TWO 
FRACTAL CONTACT AREA 
DISTRIBUTION 

A. Model 

Extensive studies have already established dif- 
ferent types of fractals and their properties 
0, 0, E3 • The statistics of overlaps between 
two such fractals is however not studied much 
yet, though their knowledge is often required 
in various physical contexts. For example, 
it has been claimed that since the fractured 
surfaces have got well-characterized self-affine 
properties, the distribution of the elastic en- 
ergies released during the slips between two 
fractal surfaces (earthquake events) may fol- 
low the overlap distribution of two self-similar 
fractal surfaces @, Q, 0]- Also, using renor- 
malisation group technique Chakrabarti and 
Stinchcombe Q have analytically shown that 
for regular fractal overlap (Cantor sets and 
carpets) the contact area distribution follows 
power law. 

Here, we study the distribution P(m) of 
contact area m between two well-characterized 
fractals having the same fractal dimension. 
We have chosen different types of fractals: reg- 
ular or non-random Cantor sets, random Can- 
tor sets (in one dimension), regular and ran- 
dom gaskets on square lattice and percolating 
clusters 0, 0] embedded in two dimensions. 
We find a universal scaling behavior of the dis- 
tribution: 

P'(m') = L a P(m,L);m' = mL~ a , (1) 

where L denotes the finite size of the fractal 
and the exponent a = 2{df — d); df being the 
mass dimension of the fractal and d is the em- 
bedding dimension. Also the overlap distribu- 
tion P(m), and hence the scaled distribution 
P'(m'), seen to decay with m or m' following 
a power law (with exponent value equal to the 
embedding dimension of the fractals) for both 
regular and random Cantor sets and gaskets: 

P(m) ~ m~' 3 ;/3 = d. (2) 

For the percolating clusters, however, the over- 
lap distribution takes a Gaussian form. It may 



be noted that the normalisation (to unity) re- 
striction on both P and to ensures the same 
scaling exponent a for both. 
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B. Overlaps between regular fractals 



Here we construct three types of regular 
fractals: regular Cantor sets of dimension 
In 2/ In 3, regular Cantor sets of dimension 
In 4/ In 5 and regular gaskets of dimension 
In 3/ In 2 on a square lattice. These regular 
fractals are constructed following a repetitive 
procedure in successive generations (n) such 
that the self-similarity is strictly maintained 
at every stage. For example, a regular Can- 
tor set of dimension In 2/ In 3 is formed in the 
n — > oo limit of the set obtained by removing 
the middle one-third portion of each occupied 
set at every generation, starting from a com- 
pact set of size L = 3™ at n = 1 (see Fig. 1). 
For the set with dimension In 4/ In 5, one sim- 
ilarly removes the middle one-fifth portion at 
each generation (here L — 5™ for finite n-th 
generation) . 

(a) 



n=l 



n=2 
n=3 

(b) 



n=3 



(d) 



Fig 1: (a) A regular Cantor set of dimension In 2/ In 3; 
only three finite generations are shown, (b) The over- 
lap of two identical (regular) Cantor sets, at n = 3, 
when one slips over other; the overlap sets are indi- 
cated within the vertical lines, where periodic bound- 
ary condition has been used, (c) A regular gasket of 
dimension In 3/ In 2 at the 7th generation, (d) The 
overlap of two identical regular gaskets at same gener- 
ations (n = 7) is shown as one is translated over the 
other; periodic boundary condition has been used for 
the translated gasket. 

To study the overlap between two such iden- 
tical fractals, the boundary effects are avoided 
using periodic boundary condition for one and 
we consider the other set to slip over the first. 
At each step of such slip or translation, we 
count the overlapping filled sites (black parts 
or dots) present in both the fractals and the 
total number of such sites gives the size of 
overlap to. Thus, if one sequentially translates 
one fractal over another, various overlap values 
(to) are obtained which in turn give the distri- 
bution P(m). Since no randomness is involved 
here, we do not need any configurational av- 
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eraging. These results are shown in Fig. 2, 
where the Cantor sets and gaskets are gener- 
ated for finite generations n: L — 3™ for Can- 
tor sets with df = In 2/ In 3, L — 5" for Cantor 
sets with df = In 4/ In 5 and L — 2™ for gaskets 
with df = In 3/ In 2. The overlap distributions 
P(m, L) are fitted to the scaling forms (1) and 
(2). The results indicate their validity in the 
large n (or L) limit. 
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Fig 2: The scaled distribution plot of P'(m') = 
P(m, L)L a against the scaled overlap m! = mL~ a 
for two identical regular fractals: (a) Cantor sets with 
df = In 2/ In 3 at various finite generations: n = 10 
(square), n = 11 (plus), n = 12 (cross) and n = 13 
(star); (b) Cantor sets with df = In 4/ In 5 for fi- 
nite generations: n = 6 (square), n = 7 (plus), 
n = 8 (cross) and n = 9 (star); (c) gaskets with 
df = In3/ln2 for finite generations: n = 7 (plus), 
n = 8 (cross), n = 9 (star) and n = 10 (square). Note 
that in all the three cases a = 2(df —d) and the dotted 
lines indicate the best fit curves of the form a(x — b)~ d ; 
where d is the embedding dimension [d = 1 for (a) and 
(b) and d = 2 for (c)]. Insets show the unsealed distri- 
butions P(m) for overlap m in different cases. 
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C. Overlaps between random fractals 



Here we construct three types of random 
fractals: random Cantor sets of dimension 
In 2/ In 3, random Cantor sets of dimension 
In 4/ In 5 and random gaskets of dimension 
In 3/ In 2 on square lattice. We can construct 
a random Cantor set of dimension In 2/ In 3 by 
removing randomly any of the one-third por- 
tion at every stage or generation (see Fig 3). 
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fractal dimension) in two different realisations. The 
overlap sets are indicated within the vertical bars, (c) 
A random realisation of a gasket of dimension In 3/ In 2 
fl=l at 7th generation, (d) The overlap of two random gas- 

kets of same dimension and of same generation but 
generated in different realisations. 

n=2 
n=3 



Here the structure of the sets change with 

* * configurations as randomness is involved. We 

therefore take the overlap between any two 

such sets at finite generation n having same 
n _3 dimension but of different configurations. We 

count the portions present in both the sets, 
and this gives the total size of overlap m. 
Clearly the overlap sizes (m) change from one 
(c) pair of configurations to other- and we deter- 

mine their- average distr ibution P(m). The fi- 
nite size L of the fractals is similarly related 
to generation number n as discussed in the 
previorrs section. The average distributions 
P(m. L) for- finite size L of the fractals are 
..J determined using 500 such configurations for 

■3 each of the three types of random fr actals. 
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Fig (a) A random Cantor set of dimension In 2/ In 3; 
only three finite general ions are shown, (b) Overlap of 
two random Cantor sets (at n = 3; having the same 
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D. Overlaps of percolating clusters on 
square lattice 



Here we study the overlap distribution of two 
well-characterised random fractals; namely the 
percolating fractals Efficient algorithms 

are available to generate such fractals. It 
seems, although many detailed features of the 
clusters will change with the changes in (par- 
ent) fractals, the subtle features of the overlap 
distribution function remains unchanged. 
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Fig 4: The scaled distribution plot of P'(m') 
P(m, L)L a against the scaled overlap m! = mL~ 
for random fractals: (a) Cantor sets with dimension 
df 



In 2/ In 3 for finite generations: n = 14 (square), 
n = 13 (star), n = 12 (cross) and n = 11 (plus); (b) 
Cantor sets with df = In 4/ In 5 for some finite gen- 
erations: n = 10 (plus), n = 9 (cross), n = 8 (star) 
and n = 7 (square); (c) gaskets with df = In3/ln2 
for finite generations: n = 11 (square), n = 10 (star), 
n = 9 (cross) and n = 8 (plus). For all the three cases 
a = 2(df — d) and the dotted lines indicate the best fit 
curves of the form a(x — b)~ d ; where d is the embed- 
ding dimension of the fractals [d = 1 for (a) and (b) 
and d = 2 for (c)]. Insets show the unsealed distribu- 
tions P(m) for overlap m. 



We generate numerically several site per- 
colating clusters at the percolation threshold 
(p c = 0.5927 12]) on square lattices of lin- 
ear size L using Hoshen-Kopelman algorithm 
[13 03| For the overlap of any two clusters, we 
count the number of sites M which are occu- 
pied in both the clusters (see Fig 5). This gives 
the overlap size m = M/L d between the frac- 
tals. As the realisations change (keeping the 
fractal dimension df of the percolating cluster 
same) m varies and we find out its distribution 
P(m,L). We find that the distribution shifts 
continuously as L increases and has a finite 
width which diminishes but very slowly with 
L. This shows, the emerging length scale asso- 
ciated with m is no longer a constant, rather 
it depends on L. This arises due to the frac- 
tal nature of the original clusters, where the 
occupation of the sites are no longer random 
events, but are correlated [T3 |. Hence, for a 
system of size L, the probability of occupa- 
tion grows as L df for any of the fractals and 
as [L d f/L d ] 2 = L 2{d s~ d ^ for the overlap set. 
If this is the origin of the L dependence in 
P(m, L), then the distributions for different L 
should follow a scaling behavior as indicated 
by eqn. (1). 
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a connected one. 

From the overlap between all the pairs of 
cluster configurations (typically around 500 
for L = 400), we determine the distribution 
P(m,L). The data are binned to facilitate 
storage and to make the distribution smooth 
(Fig 6 (a)). We have also studied the nature 
of the distribution P(m) for percolation clus- 
ters generated with lattice occupation proba- 
bility p above p c of the square lattice. Here 
the distributions become delta functions and 
the height of the delta function increases with 
system size L (Fig 6 (b)). 
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Fig 5: The overlap between two percolating clusters; 
(a) and (b) are two typical realisations of the same 
percolating fractal on square lattice (df ~ 1.89) and 
(c) their overlap set. Note, the overlap set need not be 
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Fig 6: (a) The scaled distribution plot of P'(m') = 
P(m, L)L a against the scaled overlap ml = mL~ a for 
percolating clusters grown with probability p = p c = 
0.5927 on square lattice (df ~ 1.89) for finite sizes: 
L = 100 (cross), 200 (star) and 400 (square). Here 
a = 2(df — d). The dotted line indicates the best fit 
curve of the form aexp(— (x — fe) 2 /c); where a (= 
4.5) , b (= 0.38) and c (= 0.018) are constants. Inset 
shows the unsealed distribution P(m) against overlap 
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size m. (b) The distribution P(m) against overlap size 
m for percolation clusters generated with occupation 
probability p above p c (p = 0.7 and p = 0.9) on square 
lattice [for finite sizes: L = 100 (cross), 200 (star) and 
400 (square)]. Clearly the distributions P{m) are delta 
functions and the height of the delta function increases 
with the system size L. 



III. CONCLUSION 

In the context of studying the earthquake mag- 
nitude distribution, we have studied here a 
model where the earthquake fault and the 
moving tectonic plates are assumed to be self- 
similar fractals. Essentially, we study the con- 
tact area distribution of two nominally iden- 
tical fractals. Our study on different sets of 
fractal surface overlap shows that the contact 
area distributions are not always power laws. 
Although for Cantor sets and gaskets hav- 
ing different fractal dimensions (both for reg- 
ular and random cases), the contact area dis- 
tributions have got power laws, for percolat- 
ing clusters (at percolation threshold) it takes 



a robust Gaussian form which we can't ex- 
plain even qualitatively. In particular, we find 
that in all the cases the distributions P(m, L) 
show an universal finite size (L) scaling be- 
havior P'(m') = L a P(m,L); ml = mL- a , 
where a — 2(df — d). The P{m), and con- 
sequently the scaled distribution P'(m'), have 
got a power law (2) decay with m (with de- 
cay exponent equal to d) for both regular and 
random Cantor sets and also for gaskets. For 
percolation clusters, P(m) (and hence P'(m')) 
have a Gaussian variation with m (Fig. 6a). 
It may be mentioned that, since in many cases 
the bulk thermodynamics may be mapped to 
some reduced dimensional interface problems 
and since their overlaps corresponds to phys- 
ical quantities (e.g. , the replica overlap in 
the physics of glasses ^3]), such fractal over- 
lap distributions are of quite general interest. 
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